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Abstract 



This paper is devoted to the convergence and stability analysis of a 
class of nonlinear subdivision schemes and associated multi-resolution 
transforms. These schemes are defined as a perturbation of a lin- 
ear subdivision scheme. Assuming a contractivity property, stability 
and convergence are derived. These results are then applied to vari- 
ous schemes such as uncentered interpolatory linear scheme, WENO 



scheme [13jj , Power-P scheme [16| and a non linear scheme using local 



spherical coordinates [18|. 
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1 Introduction 



Multi-resolution representations of discrete data are useful tools in several 
areas of application as image compression or adaptive methods for partial 
differential equations. In these applications, the ability of these representa- 
tions to approximate the input data with high accuracy using a very small 
set of coefficients is a central property. Moreover, the stability of these rep- 
resentations in presence of perturbations (generated by compression or due 
to approximations) is a key point. 

In the last decade, several attempts to improve the property of classical 
linear multi-resolutions have lead to nonlinear multi-resolutions. In many 
cases, this nonlinear nature hinders the proofs of convergence and stability. 

In [H , in the context of image compression, a new multi-resolution trans- 
form has been presented. This multi-resolution is based on an univariate 



nonlinear multi-resolution called PPH multi-resolution (see [121 ] in the con- 
text of convexity preserving). It has been analyzed in terms of convergence 
and stability of an associated subdivision scheme following an approach for 
data dependent multi-resolutions introduced in 0|. Due to nonlinearity, the 
stability of the PPH multi-resolution is not a consequence of the convergence 
of the associated subdivision scheme. It has been established in [2], present- 
ing the PPH subdivision scheme as some perturbation of a linear scheme 
following @, Q, and 



The aim of the present paper is to generalize the results presented in 
for a general family of nonlinear multi-resolution schemes associated to an 
interpolatory subdivision scheme Snl : Z°°(R) — * Z°°(R) of the form: 

V/ G r(M), VnGZ f S NL (f) 2n+1 =S(f) 2n+1 + F(5f) 2n+1 , 

[ ^NL(J)2n=Jn, 

where F is a nonlinear operator defined on Z°°(R), 5 is a linear and continuous 
operator on Z°°(R) and S is a linear and convergent subdivision scheme. 
Considering two subdivision schemes Snl and S, it is always possible to 
introduce the difference F = Snl — S. If one assume some properties of 
polynomial reproduction (see section [3]), as shown in F is in fact a 
function of differences, i.e. of df defined by df n = f n +i — fn- 



Theorems [T] and that are the main results of this paper, establish that 
if F, S and 5 satisfy some natural properties, then the subdivision scheme is 
convergent and the multi-resolution is stable. 

The paper is organized as follows. In section [2] we recall briefly the 
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Harten's interpolatory multi-resolution framework which is the natural set- 
ting for our work. We precise the class of schemes under consideration and 
we establish the main results in section [3l Various applications are presented 
in section SJ 

2 Harten's framework and basic definitions 

In Harten's interpolatory multi-resolution, one considers a set of nested bi- 
infinite regular grids: 

X j = WJneZ, x{ = n2~ j , 
where j stands for a scale parameter and n controls the position. 

The point-value discretization operators (or sampling operators) are de- 
fined by 

Vj ■ f e C(R) ~ f = U 3 n)n& ■= (f(xi)Uz e V 3 , (2) 

where V 3 is the space of real sequences and C (R) the set of continuous func- 
tions on R. 



The reconstruction operators IZj associated to this discretization are any 
right inverses of T>j on V 3 , that is, any operators IZj satisfying: 

(K j f)(xi) = fi = f(xi). (3) 

For any j, the operator defined by VjlZj + \ acts between a fine scale j + 1 
and a coarser scale j. Here, it is a sub-sampling operator from Vj + \ to Vj. 

The operator defined by Vj + \R,j acts between a coarse scale j and a finer 
scale j + 1 and is called a prediction operator. A prediction operator can be 
considered as a subdivision scheme [15( from Vj to Vj+i- We say that the 
subdivision scheme S defined by (f J ) t— > S(f 3 ) = Vj + i1Zj(f 3 ) is uniformly 
convergent if: 

V/ e 3f°° G C°(R) such that lim sup|^'(/) n - / 00 (2 _i n)| = 0. 

We note f°° = S°°f. 

Since for most function /, Vj + \R,jf 3 ^ f 3 , details, called d 3 and 
defined by d 3 = f 3+1 — T>j + \Rjf- 7 , should be added to Vj^TZjf 3 to recover 
f 3+1 from f 3 . The multi-resolution decomposition (see [3j, [10(], 11] for 
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precisions) of f L is the sequence {/°, d°, . . . , d L Moreover, the multi- 
resolution transform is said to be stable if: 

3C such that V/ L ,/ L ,j<L 

ll/ J '-/ J 'lloo ^cfwf-plU + £||d fc -d*||ooV) 

V k=0 J 

||/°-/°||oo <C||/ J '-/ J '||oo, (5) 

H^-^Hoo <C\\f -fW^ Vfc, 0<k<j-l, (6) 

where {/°, d°, . . . , d 1 ^ 1 } is the multi-resolution decomposition of f L . 

When the prediction operator Dj + iTZjf J is linear, the convergence of the 
associated subdivision scheme implies the stability of the multi-resolution 
analysis. In the non linear case, it is not the case and there is no general 
result for the multi-resolution analysis stability 



3 A Class of Nonlinear Subdivision Schemes 

Introducing S a linear, reproducing polynomials 1 up to degree P and conver- 
gent interpolatory subdivision scheme we consider nonlinear interpolatory 
subdivision schemes that write 

SNL(f j hn+l=S(P)2n+l+F(5P )2n+l /„n 
S N L(f j hn=fi 1 ' 

where F is a nonlinear operator defined on Z°°(Z) and 5 is a continuous lin- 
ear operator on Z°°(Z). 



3.1 Convergence analysis 

We have the following theorem related to the convergence of the nonlinear 
subdivision scheme Snl'-' 

Theorem 1 If F,S and 5 verify: 



3M>0 such that \/d G Z°° ||F(cf)||oc < M||cZ||oo (8) 
3c <1 such that \\SS(f) + 5F(6f)\ U < c\ \Sf\ U, (9) 

The interpolatory subdivision scheme S reproduces polynomials of degree P if, for 
any polynomials V of degree less or equal to P, if /„ = V(x J n ) then S(f)2n+i = ^(^n+i) 
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then the subdivision scheme Snl is uniformly convergent. Moreover, if S 
is C a ~ convergent (i.e for all f G i°°(Z), S°°(f) G C a ~ 2 ) then, for all se- 
quence f G Z°°(Z), Sfi L (f) G CP~ with [3 = min (a, log 2 (c)). 

Proof Using hypotheses (JSj) and ([9]) and the definition of Snl, we get : 

\S N L(f j )2n + l-S(f) 2n+1 \ < MW5PIU 

[|fe(/ J ')-W)l|oo < MlltfOW*- 1 )^, 
||^l(/ J )-^(/ J )IIoo < Mcll^lU, 

that can be rewritten as : 

\\S N L(f j )-S(ft)\\oo < McipfU. 

Writing : 

\\S N L(f J )-S(f J )\\oo < M||,5/ || oo 2^( c ) (10) 

the convergence of the subdivision scheme Snl can be obtained applying 
theorem 3.3 of [7[. 

In our context, this theorem applies as follows : 

If S is a linear C a convergent subdivision scheme reproducing polyno- 
mials up to degree P and if Snl is a perturbation of S in the sense that, 
calling f k := S NL (f°) for all /° G 



\\S NL (f k )- 5(^)1100 = 0(2" 



then SVl is C^ 3 convergent with /? => min(P, sl,v). 

It follows that if S is C a convergent then Snl is at least C@ convergent 
with = min (a, /002(c) )■ 

□ 



Remark 1 When F is linear, theorem [l\is a consequence of theorem 6. 2 in 

a 

2 For < a < 1, C a = {/ continuous, bounded and verifying Va < a, 3C > 
0, Vx,y £ R, |/(a:) - /(y)| < C7|as - y| a }. For a > 1 with a = p + r > 0, p e IN 
and < r < 1, C* a ~ (R) = {/ with f (p) G C r ~ } 
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Remark 2 In many of our examples, S is the two point centered linear 

scheme defined by S(fi)2n+i = 2 " +1 which limit function is in C l ~ . 
Therefore, as soon as the non linear scheme Snl verifies hypothesis and 
(0|) with c>\, the S NL is C^ l ° 92 ^~ convergent. 

Remark 3 Hypothesis can be weakened as: 

3 P £N 3c <1 such that \\8(S p NL f)\\ OQ < cWSfW^. 
The proof remains the same except that : 

\\S NL {f 3 )-S{P)\U < Mc\\8PU 

becomes : 

\\S N L(P)-SU j )\\oo < M\\5(S p NL f-v 
\\S NL (f j ) - S(f)\\oo < Mc\\8f-P 
that can be rewritten, for j = i[p], as: 

\\s NL {f) - sinu < Mc^wsfu. 

The conclusion is reached applying theorem 3.3 offji]. 

Remark 4 A straightforward generalization of theorem U\ can be obtained 
introducing two linear operator 8\ , 82 and a perturbation of the form F(8±f, 82 f) ■ 
Under the following hypotheses: 

3M > such that \F(d,d')\ < Mmax (\\d\loo, Hd'Hoo), (11) 
3c>l such that IMSuLifMoc < c max (| I loo, I l^/l |oo), (12) 
\\82(S NL (f))\\oo < cmaxdl^/IUJI^/lloo), (13) 

for all d,d' E l°°,f S l°° , the scheme Snl is uniformly convergent. 
Remark 5 We can also apply theorem [7] to bi-variate schemes written as 

S NLl (x j ,y j ) 2n+ x \ I a^+ +1 \ _ f S(x j ) 2n +i + Fi(8x j ,8y j 



S N L(x>,y>) y S NL2 (xi,yi) 2n+1 J ■ ^ ^ J y S&) 2 n+1 + F, .(,»' . 8 { f' ) 
If the following conditions are satisfied for i = 1, 2 

3M > such that \Fi(d,d')\ < M max (\\d\loo , \\d f \\oo) , 
3c> 1 such that \ \8(S NLl (x,y))\ U < cmax (| \8x\ 1^, | \8y\ |oo), 
for all d,d',x,y S l°° , the scheme Snl is uniformly convergent. 
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3.2 Stability analysis 

We now consider the multi-resolution analysis associated to the subdivision 
scheme ([7]) recalling that, for any sequence f J , the details d? are defined by 

di = Hull ~ SNL(Phn+l- 

We have the following theorem concerning the stability of the multi- 
resolution : 



Theorem 2 If F, S and 5 verify: 3M > 0, c < 1 such that V/, g, di,d%, 

||F(di) - F(da)||oo < Af ||di - dalloo, (14) 
||*0W - W)[|oo < 45(f - 5 )|U, (15) 

then the multi-resolution transform associated to the non linear subdivision 
scheme Snl is stable. 



Proof 

We first prove (jU) : 

Due to the interpolatory property we only consider |/| n+1 — /^n+il- 

Since S is a convergent linear scheme, we have, using the stability of the 
linear scheme S: 3C > such that 



l/L+l-ifn+ll < ^^ll/ -/ ||oo + Ell/ fc - 5 (/ fc " 1 )-/ fe + 5 (/ fe " 1 )llo°J 

< C ( ||/° - /°||oc + ll^" 1 + F(5f k ~ l ) - ~d k ~ x - FiSp- 1 ) 



From (0 : 



k=l 



( i-i 3 

l/L+l-^n+ll < C lir-^lloo + Ell^-^lloo + ^Ell^ 1 )-^^ 

V A:=0 fc=l 

Concentrating on the last right hand side term we get: 
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? °ioo 



fc=l 



ElW/*- 1 )-^/*" 1 )!!*, ^ \\5{f)-5{~f 

j 

+ E (\\ S ( S NLf k - 2 ) ~ ^OW^IIoc + H^- 2 - ^- 2 ||oo) 



From (11511 we get: 



k=2 



j 3-2 

Eh^ 1 ) - < y (/ fc ~ 1 )ii°c < n*(/°) - *(/°)iioc + E (W) - *(/ fe )H~ + n w * - w ~ 

fe=l fc=0 

J-2 / fc \ 

< E cfc n^° - 5 f°\\°° + E cfc ~'n^ - 5 *\\°° • 

fc=0 V 1=0 J 

Since < c < 1 we get finally: 

i-i 



ll^'-flloo < C'llf-^Hoo + C'Ell^-^ 

fe=0 



+A/C— - <S(/»)|U + £ ||id* - w 



|00 I 5 



fc=0 



and, using the continuity of 6, we get (jU) with a constant 

C = C7 , + MC'|| ( 5|| 00 



1 - c 



We now establish (JSJ) et ©. 

Equation ([5]) is a direct consequence of the interpolatory properties. 
For ©, we have, for < A; < j — 1 : 



\d k -d k \ < - f*+i - S (f k ) - Stf")]^ + \\F(5f k )-F(Sf k )\\ o0 . 



Using the property ([6]) for the multi-resolution associated to S, hypothesis 
(fTil) and the continuity of 5, we have : 



\d k -d k \ < C / ||^-^|| 00 + M||5|| 00 ||/ fc -/ fe || 00 . 



'fen 

loo 



s 



From ([5|) for the multi-resolution associated to Snl we have 
\d k n -d k n \ < C\\f- fill,, + M\\5\\ o0 \\f3- 1 -f3- 

and therefore we get © with C = C + M||<5||oo. 
□ 



Remark 6 As previously, we can again consider a weaker formulation for 
hypothesis $15\) such as: 

3peN, 3c <1 such that \\6(S p NL f - S^)||oo < c\\5(f - g)^. 

Under this hypothesis, the stability of the subdivision scheme can still be 
established. However, the multi-resolution stability is not ensured. To get it, 
a stronger hypothesis like: 

3p S N, 3c < 1, such that 

P-2 

IW/ P /-/)IU<c||<5(/- ff )|U + M^\\d k (f)-d k (g)\\oo, 

k=0 

is required. 



4 Applications 

This sections is devoted to applications of the previous results to three spe- 
cific subdivision schemes (linear and nonlinear) available in the literature. 
We provide for each of them, the proofs of convergence and stability. 
In all this section, given f 3 = {f 3 k )k^2Z we note: 

df = (C)ne^ with dfi = fi +1 -fi, (16) 
Df = (DfDnezz with Dfl = f n+1 -2fl + fi_ 1: (17) 

(18) 
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and, more generally D l f j = (D l fi) n( zzz with: 



D l fi = D(D l - 1 fi) n = J2(-l) i Cyn-i+i with Cl = — — - .(19) 

i=0 ^ 

4.1 Multi-resolution analysis associated to a linear fully non 
centered Lagrange interpolatory subdivision scheme 

As it has been said before, for linear scheme, the stability of the multi- 
resolution analysis is a consequence of the convergence of the subdivision 



scheme (see [111 ] ). Therefore, we only consider here the convergence of the 
subdivision scheme. 

The convergence of centered linear interpolatory schemes is well know 
since Delauriers and Debuc [6fl. 

For linear but non centered schemes there is no general results of conver- 



gence. Moreover the general tools proposed in 15)] are very fastidious to 
apply and don't provide general results. 

In this subsection, we focus on completely decentred Lagrange interpo- 
latory linear schemes. In order to apply our theoretical results, we consider 
S the two point centered linear scheme and express any right hand side 
excentred scheme Sp (where P stands for the number of point of the con- 
sidered stencil) as a perturbation of it. Precisely, if we write Sp{f 3 ) 2n+ \ = 
S(f j ) 2n +i + Fp(5pf j ) 2n +i we get: 

If P is even, 



P-2 

F P (5 P P) 2n+1 = + ^fn+l+l 2^{k-mk+l)\ 

k=2 keven 
P-3 

,^±i f (4fc+5)(2fc-l)! 



k=l kodd 



2 2fc+l(fc_l)!(/ c+2 )! • 



and,if P is odd 



P-3 ^ 

F P (5pP) 2 n+l = + Yl 2 a*(fc-l)! 1 (fc+l)! 

k=2 k even 
P-4 



E n i£±I , (4fc+5)(2fc-l)! 
u J n+^ 2 2fc + 1 (fc-l)!(fc+2)! 



k=l kodd 
P-l 



r>^ f (2AT-3)! 
u J n+*f± 2 2(N-2) (jv-3)!(AT-1)!- 
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Numerical evaluation of the perturbation terms for 4 < P < 9 is given in 
tabled) 



F P (Df) 2 



n+1 



j§Df n+ i + ^D/ n+2 

■J§Df n+ i + j^Df n+2 — jfg-D 2 /rt+2 
16-D/n+l + J%Df n+ 2 — ^D 2 f n+ 2 + ^fg-D 2 / 

-3 7-i i , 1 t-> x 17 n2f , 7 n2 . 



10 JI 'Ti ' 10 J "T^ ZOb J "T^ ' ZOO J 

J§Df n+ i + jjrDfn+2 — 2lk-D 2 fn+2 + 2EE ^ f n + 3 ~ JM^D 3 f n+3 

if -D /n+l + yg£>./n+2 - 2lk D2 f n + 2 + f n +^ ~ 204$^ -^+ 3 + 2Si8 D ^ fn+i 

if -D /n+l + jQ D f n+ 2 ~ ^D 2 f n+2 + 25g J D 2 / n+ 3 - -gTjfg D 3 /„ +3 + 2518 -^A^ 
_ 429 r)4 <• 



429 n 4 f 
"32768-^ ■' n + 4 



Table 1: Perturbation term Fp(Df)2 n +i for different values of P 



It also appears that Sp can be written naturally as a perturbation of Sp_2, 
for even values of P and as a perturbation of Sp-i for odd values of P. 
Indeed, we have: 
When P is even: 



Sp{f)2n+1 - Sp- 2 (f)2n+l + 2 2(P-2)(p_3 ) ) |'(p_ 1 )i D 2 f n+ P 

(4P-8)(2P-7)! n^f p 

2 2F-5(p_4)!(p_ 1 )! ^ 7 n +£_l) 

(20) 

and when P is odd: 

Sp(f)2n+ 1 = 5 , p_i(/) 2 n+l - 2 2(P-2)(p„3)i'(p_ 1 )i D ~ f n+ E=±- 

2 (21) 

In both cases, its is easy to check that the function F defined by F = 

P — 2 P — l 

F(D~z~ f) when P is even and by F = F(D~z~ /) when P is odd, is linear 
and continuous. Therefore, the convergence can be reached as soon as the 

P — 2 P—l 

contractivity hypothesis © for or is satisfied. 

Direct calculations provide the estimates gathered in table (2) It then fol- 
lows from theorem [T] that all the fully excentred interpolatory subdivision 
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p 


perturbation term 


contractivity estimate 


4 


F(Df) = -±Df n+1 + ±Df n+2 


\\D(Sif)\\oo < IIP/Hoc 


5 


F ( D2 f) = -lk D2 fn + 2 


||£> 2 (5 5 /)||oo<ilP 2 /l|oo 


6 


F(D 2 f) = -j^D 2 f n+2 + 25QD 2 f n+3 


|P 2 (5 6 /)|U < ^H^/lloo 


7 


F (D 3 f) = ~^D 3 f n+3 




8 


F(D 3 f) = _ 2048- D3 /n+3 + 2m DZ fn+i 


IP 8 (5 8 /)||oo<SIWIIoo. 


9 


F(P f) ~ 32768 D f n + 4 


iiD 4 (5 9 /)iioo<ifii^ 4 /iioo. 



Table 2: Perturbation term (see 1201 and I20p and contractivity estimate for 
different values of P 



schemes for P < 8 points converge. 

The following comments can be made for the other situations: 

Figure [Q represents the completely decentred 9 and 10 points scheme 
iterated functions at scale 8 starting from / . From the zooming in the 
oscillating region, one can guess that the 9 point scheme converges while the 
10 point doesn't. In fact, following [l5| . and using the so called iterative 
formalism one observes numerically that the spectral radius of the iterated 
matrix for 10 points overshoots the critical value 1 while the spectral radius 
of the iterated matrix for 9 points doesn't, that confirms the guess. 

Obviously, for these linear schemes, theorem [2] applies as soon as theorem 
[Udoes. Therefore stability is ensured for P < 8. 



4.2 The 6 points WENO subdivision scheme 

WENO subdivision schemes [131 ]. are constructed using convex combination 
of different interpolatory polynomials of fixed degree. For degree 3 and 
therefore a 6 point stencil, the WENO — 6 subdivision is given by: 

dweno\J )2n+l — ^g/„_2 ^g Jn4 t l i + g JJn 

5a + 2a x v f j 5a + a x } oq j 

+ V 1 + g )Jn+l Jn+2 "r J n+3 
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0.5 

(a) 9 point scheme 




0.45 0.47 0.5 

(c) 9 point scheme (zooming) 



3 
2 
1 


-1 
-2 
-3 



3 

2 

1 


-1 
-2 
-3 



0.5 

(b) 10 point scheme 




0.45 0.47 0.5 

(d) 10 point scheme (zooming) 



Figure 1: Iterated function at scale 13 starting from / 8 denned by /| = 1 if 
k2~ s > .5 and /? = elsewhere 



where the coefficients on control the convex combination and therefore sat- 
isfy oti > and ao + «i + «2 = 1- 



In [13J, these coefficients are defined as: 



with 



ao + «l + 0.2 



d; 



1 (e + 6,) 2 

where bi, defined as a function of the first difference df is an indicator of 
smoothness while dj and e are fixed positive constants. A set of possible 
values for these constants is suggested in [171 ]. 
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The convergence of the associated subdivision scheme has been stud- 
ied in 0] . We present an alternative proof, using theorem [TJ 

First, the WENO — 6 subdivision scheme is written as a perturbation 
of the linear two point interpolation scheme as: 

Sweno(f j )2n+l = ~ + fl+2 ( 22 ) 

3ao j-oi j Qi + 3a 2 ■ a 2 j 
16 16 Jn iq Jn-l 

with «o = o:o(df 3 ), Oil = a\{df 3 ) and a 2 = a2(df J ). 
We then have the following proposition : 

Proposition 1 The WENO — 6 subdivision scheme is convergent and, for 
any initial sequence f 3 , the limit function belongs to C~ l ° 92 ^ . 

Proof 

According to remark [U the proof can be performed in three steps con- 
sidering that F is a function of df 3 and Df 3 . 

First, according to the definition of F and to the properties of a^, we 
have 

\F(d,D)\ < imaxdldlUHDlDoo). 

Second, we prove (jlip for the first difference operator d: 
We have, for f € l°° : 

d(S weno (f))k — 5'uieno(/)fe+l <S'toeno(/)A; 

We have to consider two cases, according to the parity of k. We give the 
details for k = 2n + 1, the even case being similar. 
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d{S W eno{f))2n+l — S weno (f)2 n +2 ~ S weno (f)2n+l 



, fn + fn+i a n f 3a +ai nf 

Jn+l 7, -T^ 1J Jn+2 H 77 ^JnH 

2 lb lb 

lb lb 

dfn ao n , , 3a + «i 

— - -D/ n+2 + 16 Df n+1 

Ql + 3a 2 a 2 

H 77 ^ n ~~ T^ v J n - 1 

lb lb 



Since ao + ai + a 2 = 1 and < «i < 1, we have: 

KW/)Wi| < |||4T||oo+ (S + ^ + S + S + tI + S) 1 ™ 1 -' 

< P4f||00 + —^110/1100, 

< llM/lloo + lll^/Noo, 

< Imaxdld/lloclp/Hoo.) (23) 

Third, we prove inequality (|9|) for the second difference operator D. 
Again, two cases have to be considered: 
• For k=2n+l, then: 



D(S weno (f))2n+l — fn+l — 2S weno (f)2n+l + fn-i 

a 3a + ai ai + 3a 2 a 2 

= — -TT-^/n+2 H £ Ujn+l H 5 -^Jn 



Using < ao < 1, < ai < 1, < a 2 < 1 and ao + ai + a 2 = 1, we get: 
|-D(5 , 1i , eno (/)) 2 n+i| < — rr -||-D/||oo, 

o 

< ip/Hoo. (24) 
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• For k=2n, then : 



D{S weno (f))2n — Sweno\f)2n+l 2/ n + S weno (f) 



weno\J )2n—li 
fn+1 — %fn + fn-1 «0 n , 

2 + 16 /n+2 

2a + ai 3a + 2ai 3a 2 n , 

77 UJn+i — h —r^-lJJn 

16 lb 16 

16 lb 



Using a + ai + a2 = 1, we get: 



D(S weno (f))2n — TT^-D/n + 77-D/n+2 77 ~Dfn+l + —^Df n 

16 16 lb 16 

77 Ujn-l + —Uj n -2 



Then, with < ckq < 1, < a± < 1 and < a 2 < 1 : 

< (^ + 3ao+3 1 a 6 1+3a2 )i^/iu 

< 5P/II00. (25) 



Therefore, from (|23l) . (I24p and (125j) . we obtain the inequality: 

3c < 1 V/e r max(||d(^ eno (/))|| 00 ,|| J D(5 weno (/))|| 00 ) < ~ max (||d/||oo, ||Af||oo). 

Finally, using remarks [Hand [2] we get the convergence of the WENO — 6 
subdivision scheme to a C~ l ° 92 ^ function. 
□ 
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4.3 Power-P subdivision scheme: definition and convergence 



In the same vein as the PPH scheme ([H), t ne power P scheme is a 
four point scheme based on a piecewise degree 3 polynomial prediction. 
Considering Sc the centered four point Lagrange interpolation prediction 
that reads: 

(5r(/)Wi = fjl± Y^ ~ l Dfn+1 2° fn > (26) 

the definition of the Power-P subdivision scheme is based on the substi- 
tution of the arithmetic mean of second order differences, EIi±LL£li ; by a 
general mean power p (Dfj, Dfj+x) defined in [3] for any integer p>l, and 
any couple (x, y) as: 



sign(x) + sign{y) x + y , 
power p (x,y) = — 1 



x-y 



x + y 



(27) 



Note that it coincides for p = 1, with the arithmetic mean and for p = 2 
with the geometric mean. The Power-P subdivision scheme then naturally 
appears as a perturbation of the linear two point interpolation scheme since 
it is defined by 

Spowerpif^hn+i = ^ n+ ^ n+1 - lpower p (D f w D f n+1 ). (28) 

Before establishing the convergence theorem we first prove the following 
lemma: 

Lemma 1 For any (x,y),(x' ,y') £ M 2 , the function power p satisfies the 
following properties : 

1. power p (x,y) = power p (y,x) 

2. power p (x,y) =0 if xy < 

3. power p (—x,—y) = —power p (x,y) 
4- \power p (x, y)\ < max (|x|, \y\) 

5. \power p (x,y)\ < pmm(\x\,\y\) 
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Proof 



Claims of 1 — 4 are obvious; 
Inequality 5 comes from the equality 



power p (x, y) 



□ 



sign(x) + sign{y) 



min(x, y) 



1 + 



x-y 



x + y 



+ 



We then have the following proposition: 



+ 



x + y 



p-i 



Proposition 2 The Power P subdivision scheme is uniformly convergent 
and, for any initial sequence f 1 the limit function belongs to C l for p < 4 
and Cr^D" for p>5. 



Proof 

Here again, the hypotheses of the general theorem Q] must be checked: 
We first check hypothesis ([8]) . Using property 4 of lemma [H we obtain 
for del 00 : 

\F(d)\ < ^max(|d n |, |d n+ i|) 
o 

\F(d)\ < I 



Then we consider hypothesis ([9]): 
We study as before two different cases : 

• For k = 2n + 1: 

D(Sp OU)erp (/))2n+l — — fn 2S'p OU)er p(/)2n+l ~\~ fi 



n+1 



/n+i + fn- 2 - - n+ + 2-power p (Df n , Df n+1 ) 
jpower p (Df n , Df n+ i) 
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Prom property 4 of lemma [T] we get: 

1 



D(S power (f)) 2n+1 \ < -p/lloo- (29) 



• For k = In: 

D(Sp 0werp (f^2n — "Spouier-p (/)2n— 1 2_f n -|- S powerp (^f^2n+l 
fn + fn+1 1 



2 

fn-l + fn 1 



2 

Dfn 1 



power p (Df n ,Df n+ i) - 2f n 
■power p (Df n -i,Dfn) 



(power p (D f n , D f n+1 ) + power P (D f n -i,Df n )) 



For p > 5, from property 4 of lemma [T] we get: 

3 



^SW^OOM < tIIWIIoc- (30) 



4 



For p < 4, we note D(S powerp (f)) 2 n = Z(Df n ,Df n+1 ,Df n ^ 1 ) with 

Z(x,y,z) = - - -(power p (x,y) + power p (x, z)) 
From definition 1271 and property 4 and 5 of lemma [H we have, 
if x > 0, 

| - ^(max(x,y) +max(x,z)) < Z(x,y,z) <| 

— < Z(x, y.z) < — 
4 ~ v ; ~ 2 

0< \Z(x,y,z)\ < -\x\ 
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if x < 0, 

-< Z(x,y,z) < - + -(znin(\x\,\y\) + mm(\x\,\z\)) 
|< Z(x,y,z) <(|--)|x| 
0< |Z(x,y,z)| < -\x\ 

Finally, 

|£'(^er p (/))2n| < ^P/Uoo • (31) 
From ([29]), j30l) and ([31]), we obtain: 

||£-Wr P (/)lloc < ^P/IU /orp<4. 
pS P ™er P (/)||ac < IWDfWoo forp>5. 



Finally, theorem [T] and remark [5] provides the convergence to a C 1 if p < 4 

and Cr^d)" (R) if p > 5. 

□ 



4.4 The convergence of a non linear scheme using spherical 
coordinates 

The non linear subdivision scheme studied in this section is defined in 
where it is considered as a non regular interpolatory subdivision scheme 
using local spherical coordinates. Here, we consider it as a regular subdivi- 
sion scheme applied to the M 2 point sequence Pn(xn, fn) t n& 2z- The resulting 
scheme reads (see (l8|): 



j'+i 

x 2n+l 



2 . 



4 



cos f 0& + /i(a£ 
sin [9 J n + /i(a^ 



cos 



szra 



V i + M# 



n+l 
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with : 



fuY: (33) 

»i ■ arctan I J " + J ' I. (34) 



- x J n ) 


> 2 + Un+l - 


fn+1 


- fLi \ 


X n+1 


~ X n-l) 


fn+1 


-fi\ 


x n+l 


- xi J 



-f 3 n = arctan -? j , (35) 

\ x n+l X n / 

< = ii-el (36) 

Pi+i = (37) 

and, ei, 4 e [-lirll- 

As explained in [18j], the design of /i, is performed to produce regular limit 
functions. It is then defined C 1 function that is contractive for small 

values of a and that coincides with identity for large value of a. Note that 
h = provides the classical linear two point centered scheme. 

In our context, we will note this scheme S sp hericah and S±, S 2 will stand for 
the schemes associated to each coordinates : We then get 

Si(x,f) 2n +i = Xn + ^ n+1 + (Fx(dx, df)) 2n +i, 

S 2 (x,f) 2n+ i = ^ + / n+1 + (F 2 (dx,df))2n+l, 



with 

j 

(Fi(c*r,4f)Wl = -J [cos (ft + KO) - cos [0> n+1 + h(pi +1 ))) , 
(F 2 (dx,df)) 2n+1 = ^ (sin {9i + h(ai)) - sin (V n+1 + h(^ +1 ))) . 

From ([31 [Ml ESJ, r n , n and 7 n can be written using the first divided 
difference {df 3 ,dx ] ) as: 

Tn = \/ (dxi) 2 + (dfl) 2 

' dfl + dfl_A 



6 n = arctan 
7„ = arctan 



dxi + dx 3 n _ x ^ 

' dft\ 
dxi ) 



21 



as well as a n and j3 n thanks to ([36]) and ([37]) . 
We then have the following proposition: 

Proposition 3 The scheme S sp hericai defined in (fffffj) is convergent. 
Proof 

We again check the hypotheses of theorem Q] generalized to M 2 according to 
remark We have, 

max (\dx n \,\df n \), (38) 

and therefore, for i = 1,2: 

m(d X ,df)) 2n+1 \ < 2 ^^(\d-n\,\dfn\) ^ 



a/2 

< — maxdldxlloo, ||d/|| c 



that shows that the hypothesis ([8]) of theorem [T] is satisfied. 



We now check hypothesis ([9]). 
For / 6 l°° we have 

d(Si (x,f)) 2n = Si(x,f) 2n +i - Si(x,f) 2n 

_ X n + X n +1 + rn ^ cos ^ + _ CQS (g n+i + h((3 n+1 ))) - Xr , 

= X " +1 2 + -j (cos (0 n + h(a n )) - cos (0 n +i + h((3 n+ i))) , 



and therefore 



Using the definitions of a n and /3 n we get 



KSl(*,/))2n| < + VSlDaxdldxIloo, ||4f Hoc) w + _ _ + _ ^ 
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and 

|d(S' 1 (x,/)) 2n | < 
< 



I |oo , 

max (Halloo, \\df\\ 

c 



max (l — h'(x)) \0 n - n +i\ 

xG\— 7r,7rl 



1 V^TT \ 

- + — — max |1 - h'{x)\ max (| \dx\ |oo, | | c 

2 4 l£ -7T,7r / 



The contractivity hypothesis (JH) is therefore satisfied as soon as Vx G 
[— 7r, 7r], 1 — < < 1 + For instance, the following function h: 



h(x) 



if —TT < X < 



0.55x if 

X ~ T25rf ( X ~ 2 )2 ~~ 625tt 2 ( X — 2 )( X _ y) Y < X < 2 



f < X < ~ 

— < T- < — 

7 — x — 7 



if I 



< X < TT 



which is in agreement with the criteria proposed in [18| leads to a scheme 
satisfying ©. 

Since the same sketch of proof also provides the contractivity for \d((S2(x, f))2n\ 
we get the convergence applying theorem [TJ 
□ 



5 Conclusion 

We have formulated convergence and stability conditions for non linear sub- 
division schemes and associated multi-resolutions. These conditions deal 
with the difference with a suitable linear and convergent subdivision scheme. 
Many examples show that this formulation lead to simple proofs of conver- 
gence and stability. 
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